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Аннотация 

Введение. Для улучшения диагностики заболеваний коленного сустава необходимо повысить качество обработки 
рентгенографических изображений, т.е. дать специалистам более точную информацию для анализа патологии. 
Цель исследования — показать возможности нечеткой логики в совершенствовании алгоритма определения 
опорных линий и углов сгибания колена. Для этого необходима программа, которая анализирует рентгеновские 
снимки. Известные на сегодня методы, описанные в научной и прикладной литературе, недостаточно автомати- 
зированы. В ряде случаев ортопедам и хирургам приходится вручную дорабатывать изображения, корректиро- 
вать линии. Этот пробел восполняет представленная работа. Описан созданный автором алгоритм, который не 
предполагает участия человека, автоматически идентифицирует линии и углы сгибания колена. По результату, 
выданному системой, врач может, во-первых, судить о наличии патологии. Во-вторых, сведения, предоставляе- 
мые программой, позволяют точнее планировать, проводить операции и назначать терапию. 

Материалы и методы. Использовались снимки двух рентгеновских аппаратов, которые работают в больнице 
Аль-Базель (Латакия, Сирия). Для программной реализации алгоритма задействовали язык «Питон» (Рућоп). Ре- 
шение протестировали на 500 пациентах больницы Аль-Базель. Сравнивались результаты, которые сгенериро- 
вала новая система и предшествующие версии программ обработки рентгеновских снимков. 

Результаты исследования. Создан, описан и реализован на практике алгоритм построения опорных линий и 
углов для обработки рентгеновских снимков коленного сустава. Показаны возможности нечеткой логики в авто- 
матизации обнаружения двойного порога при выявлении границ кости на изображениях. Описана работа усовер- 
шенствованного гауссовского фильтра, предназначенного для обработки рентгенограмм. 

Модифицированный метод анализа рентгеновских снимков коленных костей включает разработку алгоритма для 
автоматического обнаружения структур и аномалий в коленных суставах, определения и измерения анатомиче- 
ских параметров, оценку степени повреждения и т.д. 

Метод определения границ контуров на рентгенограммах объединил детектор Кэнни, алгоритм водораздела и 
нечеткую логику. Программа реализована в медицинской практике и показывает точность 98 %, затрачивая на 
обработку снимка менее 20 секунд. 

Обсуждение и заключение. Новая система дает высокую точность, приемлемую оперативность и не требует 
ручной корректировки снимков. Специалисты получили возможность выявить малозаметные индикаторы нару- 
шений. Кроме того, новый метод позволяет разобраться в сложных случаях, когда сочетаются несколько факто- 
ров, указывающих на возможную патологию. Широкое внедрение метода повысит качество медицинских услуг 
в ортопедии. Следует продолжить научные изыскания в данном направлении для расширения набора стратегий 
лечения заболеваний опорно-двигательного аппарата. Предстоит создать решения с абсолютной точностью, бо- 
лее высокой оперативностью обработки, а также методы, подходящие для анализа других суставов. 


Ключевые слова: опорные линии в рентгенографии, углы коленных суставов в рентгенографии, алгоритм 
Кэнни, усовершенствованный фильтр Гаусса, алгоритм водораздела 
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АБѕќгасё 
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меге сотрагеаӣ. 
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рагатеегз, аѕѕеѕѕтепі оҒ ће аӢергее оЁ іатаре, ес. Тһе теоа Юг ӣеіегтіпіпо Фе сопіоџг Боџпӣагіеѕ оп гайіоргарћѕ 
сотбіпеа ће Саппу аеѓесіог, ће уғаѓегѕһеа а]еогИт, апа Ёл72у Іоріс. Тһе ргоргат Ваз Бееп ітріетепќеа іп тейіса] 
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Введение. Анализ рентгеновских изображений позволяет выявлять и оценивать различные патологические 
состояния [1]. Значительный вклад в повышение четкости и точности интерпретаций рентгеновских снимков 
внесли алгоритмы сегментации изображений и определения границ кости. Однако такая визуализация все еще 
нуждается в улучшении, т. к. бывает не вполне определенной, оставляет возможности для разных толкований [2]. 

Достижения в области интеллектуальных вычислений позволяют задействовать нечеткую логику для устра- 
нения неточностей и двусмысленности при анализе медицинских изображений. Нечеткая логика имитирует 
мышление человека, и это надежная основа для создания алгоритмов, способных эффективно обрабатывать и 
анализировать рентгеновские изображения [3]. 

В работе [4] референтные линии и их взаимосвязи изучались с целью коррекции формы ног. В [5] клинические 
результаты лечения многоосных деформаций нижних конечностей сравнивались с тем, как это отражается на 
рентгеновских снимках. В литературе описаны программные продукты, способные анализировать рентгеновские 
снимки костей. Однако все эти системы имеют общий недостаток: иногда они неверно распознают контуры. Если 
такое происходит, то специалист дорабатывает изображение вручную. Эта проблема решается в представленной 
научной работе. Предлагается новый алгоритм анализа рентгеновских изображений с использованием нечеткой 
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логики и усовершенствованных традиционных подходов. Решение исключает необходимость ручной доработки 
снимков. Программа работает с высокой точностью и достаточной оперативностью. Она позволяет точнее выяв- 
лять и характеризовать аномалии на рентгеновских изображениях, то есть улучшает диагностику. 

Материалы и методы. При выполнении данной научной работы автор исходил из того, что оценка деформа- 
ций конечностей требует анализа анатомических и механических осей длинных костей, а также углов между 
ними. На рис. 1 показаны исходные пределы углов для здоровых людей разного возраста. Если показатель вне 
этих пределов, значит, у пациента деформирована конечность и требуется хирургическое вмешательство. 


79-835 
ий 85—90° 


а 1/ 


85—90° 


а) б) 
Рис. 1. Опорные углы между механическими осями для здоровых людей: а — для молодых; 6 —ДЛЯ ПОЖИЛЫХ [4] 


В ортопедии понятия механической и анатомической оси имеют решающее значение для представления о 
расположении и функции костей, особенно в нижних конечностях. Анатомическая ось проходит по центру кост- 
ного стержня, через середину диафиза. Механическая ось отражает направление действующих на нее сил. В ниж- 
них конечностях она проходит от центра тазобедренного сустава до центра голеностопного сустава. 

Исследуя механическую ось, можно понять несущую способность и распределение нагрузки в кости. Это осо- 
бенно важно для диагностики и лечения остеоартритов, переломов и для определения положения имплантатов 
или протезов. 

Необходимо знать особенности анатомической оси, чтобы отличить естественную кривизну от патологиче- 
ской. Это важно в том числе при оперативных вмешательствах, поскольку помогает хирургам понять, правильно 
ли расположена кость. 

Для оценки анатомии бедренной и большеберцовой костей особенно важны два показателя: 

– АГОРА (англ. аха] [епо@ оЁ ће 415а] Ғетога1 ахіѕ — осевая длина дистальной оси бедра, ОДДОБ); 

– МРТА (англ. шефа] ргохипа| ап е оЁ ће быа — медиальный проксимальный угол большеберцовой 
кости, МПУБК). 

Они крайне информативны, если нужно выявить деформации колена. АГОЕА определяет ориентацию ди- 
стальной части бедра относительно ее продольной оси, МРТА показывает наклон медиальной поверхности боль- 
шеберцовой кости относительно ее механической оси (рис. 2). 


Рис. 2. Углы между бедренной и большеберцовой костями, используемые для оценки анатомии и выявления патологии: 
1 — медиальный колено-диафизарный угол (МКДУ) = 130° (124°—136°); 2 — медиальный проксимальный бедренный угол 
(МПБУ) = 84° (80°—89°); 3 — ОДДОБ = 81° (79°—83°); 4 — МПУБК = 87° (85°—90°); 
5 — латеральный дистальный большеберцовый угол (ЛДБУ) = 89° (85°—92°) [5] 
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Углы, показанные на рис. 2, в комплексе демонстрируют возможности выравнивания коленного сустава или 
его деформации. Ортопеды ориентируются на норму (цифры в скобках), чтобы сделать выводы о типе, тяжести 
деформаций и разработать план лечения [6]. 

Алгоритм распознавания контуров Кэнни — это надежный метод, широко используемый для идентификации 
контуров на изображениях. Он строит опорные линии и углы по рентгеновским изображениям следующим обра- 
зом. На первом этапе в используется улучшенный фильтр Гаусса, который сглаживает изображение, уменьшает 
шум [7]. Далее — вычисление градиентов с использованием операторов Собеля или Прюитт. Метод часто ис- 
пользуется для определения краев изображения. Это позволяет учитывать и слабые границы, и сильные. Так ал- 
горитм точнее определяет границы и эффективнее снижает уровень шума. 

Несмотря на свою вычислительную сложность, алгоритм Кэнни остается базовым для компьютерного зрения 
благодаря способности справляться с шумом и четко различать границы. 

Для верного выбора значения стандартного отклонения сравниваются абсолютные значения точки х(ї, /) в 
сером пикселе [8]. Если |х(ї, /)| < |х(7, )| х ХС, 7) – н, то стандартное отклонение с = 1. В этом случае точка менее 
подвержена воздействию шума. Стандартное отклонение с = 1,6 указывает на значительный шум. 

Для автоматического определения границ (71, 72) Кэнни использует двойную пороговую обработку, класси- 
фицирует пикселы на изображении по интенсивности серого цвета. Вычисляем среднее значение всех пиксе- 
лов (Т1), после чего ТІ используется для разделения изображения на две области — объект и фон. Значения пик- 
селов, превышающие начальное пороговое значение ТІ, принимаются за объекты, а остальные — за фон. Теперь 
есть два набора значений: одно для объекта (ор), а другое для фона (Ро). Для каждой из двух областей (объект, 
фон) рассчитывается пороговое значение. Вычислим 72. Начнем с того, что: 

ТО = (Тоь + ТЬе) / 2. (1) 


Новое значение 70 принимается за ТІ, и процесс повторяется до тех пор, пока 70 ~= ТІ. Конечное значение 70 
принимается за 72. 

На последнем этапе применяется привязка границ для определения связанных компонентов границ на изоб- 
ражении УМЗА (англ. уа16иѕ Кпее ѕһай ап2]е — вальгусная надколенно-бедренная суставная ось). Это помогает 
объединить разрозненные границы в более крупные структуры, что улучшает восприятие и анализ изображения. 
При этом целесообразно интегрировать принципы нечеткой логики и традиционные программные подходы к 
определению границ. В этом случае используются функции принадлежности, которые присваивают степени при- 
надлежности пикселам на основе значений интенсивности серого цвета, что позволяет получить детализирован- 
ное представление о краях. 

В отличие от бинарного подхода метод нечетких множеств облегчает постепенный переход в диапазоне зна- 
чений принадлежности от 0 (полное отсутствие принадлежности) до 1 (полная принадлежность), фиксируя не- 
точности границ. Фаззификация преобразует четкие данные изображения в нечеткие множества с учетом неопре- 
деленности. Дефаззификация преобразует нечеткие результаты в четкие для дальнейшего анализа. 

Заметное преимущество метода — его адаптивность к реальным сценариям, особенно с различной степенью 
шума и неоднозначности. Однако для достижения оптимальной производительности может потребоваться тща- 
тельная настройка параметров, поэтому крайне важно найти баланс между чувствительностью и специфично- 
стью в приложениях для определения границ [9]. 

Можно добиться более детального представления о краях, если устранить неоднозначность границ изображе- 
ний. Для этого необходимо присвоить пикселам степени принадлежности. 

Каждый пиксел обрабатывается с использованием нечетких правил. Если выполняется какое-либо из нечет- 
ких правил, показанных на рис. 3, то пиксел обрабатывается как ребро, в противном случае он не учитывается. 
Это позволяет выделить границы объектов, которые отличаются от шума и других элементов изображения. Важ- 
ность данного решения доказана применительно к анализу плечевой кости по рентгеновским снимкам [10]. 


4=1 и 13=1,14=1,15=1, 
1&&12&83& 8158 81688785818 = 0 1&812&816&8 18818 = 0 

5=1 п=ЬР=ЕЗ=Т, ЕЕ 
114&124:&138:8148:&168:817&&18 = 0 | 488158 81688785 818 = 0 

6=1 15 = 1,16=1,17=1, 

1&&12&813& 814858158 817 в 818 = 0 8 11:&128:&13&81448:818 = 0 ЕЕ 
8=1 П =1,17=1,18=1, 
114&124:&138:8148:&158:816&&17 = 0 Ш 118&128:813 88148818 = 0 


Рис. 3. Правила нечеткой логики 
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Алгоритм сегментации по водоразделам используется в компьютерном зрении для распознавания объектов. 
Изначально его разрабатывали для картографии. И сейчас он воспринимает изображения как топографический 
ландшафт, где интенсивность оттенков серого соответствует высоте местности. Линии водораздела очерчивают 
отдельные объекты как границы между водосборными бассейнами. Заполняя ландшафт от локальных миниму- 
мов, алгоритм разделяет объекты по интенсивности серого цвета. 

Алгоритм определения водораздела эффективен для некоторых приложений, однако он бывает чувствителен 
к шуму и может привести к чрезмерной сегментации. На практике для повышения производительности исполь- 
зуются тщательная предварительная обработка и управляемые маркерами варианты водораздела [11]. 

Результат работы алгоритма \!айегВе4 показан на рис. 4. 


# Беѓегепсе Ипез оп Х-гау ітадеѕ 
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] | | 


0 200 400 600 800 1000 


а) 6) 
Рис. 4. Результат работы алгоритма водораздела: а — исходное изображение; б — результат определения водоразделов 


Результаты исследования. После использования алгоритма водораздела рисуются контуры для исследуемых 
точек [12]. Для построения опорных линий рассчитали крайние точки по контуру метки (выбранной области) 
кости. Затем вычислили центр масс, преобразовали изображение в двоичный формат и нашли его центр. Опре- 
делив крайние точки с четырех сторон, взяли середину вершины верхней кости, середину основания нижней 
кости и провели линию. 

Используем формулу для вычисления угла между двумя прямыми и три координатные точки: аа, Бр и сс. Две 
являются крайними точками обеих прямых, а одна — общей точкой их пересечения [13]. 


= -ахх+ ул, (2) 
х= у -В-а, х2 = уз -Ь-а, (3) 
а=1,0х(у = у) +10х(0 а). (4) 


На рис. 5 показан результат построения опорных ЛИНИЙ и углов на рентгеновских снимках. 


# Кеѓегепсе Ипез оп Х-сау ипадез = х # Кеїегепсе іпеѕ са Хау иладех == 


рові ітэде нозі ипаде Кета В ~] оз |У розй итаде Иов тэзе Кете ву Һызбо | 9 
9 2 о 2 


Тпевоібе омет в Оррег 55 - геа Тћећоіё мег в _ Оре 5 - Кее 


Рот стене а м 6 м Роу об тегет 


бее воде 2 м Я Зее эта 


Саппу Ее 7 Сатлу Едзе 


1 ппаі 


Рис. 5. Результаты использования усовершенствованного алгоритма водораздела на рентгеновских изображениях: 
а — построение опорных линий; 6— расчет углов 
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Работу детектора границ Кэнни [14] сравнили с результатами, которые показал новый, предложенный авто- 
ром алгоритм. Очевидна более высокая производительность последнего решения. Это объясняется, во-первых, 
улучшением оператора Кэнни. Во-вторых, новый метод интегрирован с усовершенствованным фильтром Гаусса, 
поэтому отпадает необходимость в дополнительном фильтре и сокращается время выполнения алгоритма. Его 
реализация значительно повысила точность и дала лучшие результаты, особенно для изображений радиальных 
линий с высоким уровнем шума. 

Новое решение внедрили в больнице Аль-Базель в Латакии. Использовались изображения с двух рентгенов- 
ских аппаратов. Программу применили к рентгеновским изображениям 500 пациентов. В 490 случаях не потре- 
бовалась ручная доработка, т. е. точность результатов — 98 %. Обработка изображения занимает около 20 секунд. 

Предложенный алгоритм программно реализован на языке «Питон» (Руѓћоп) с графическим интерфейсом. 
Схема создания опорных линий и углов для обработки рентгеновских изображений коленного сустава призвана 
улучшить интерпретацию рентгенографических изображений (рис. 6). 


Начало 


Исходное изображение 


Алгоритм Кэнни 
(автоматическое 
вычисление ТІ, Т2) 


Морфологические 
операции 


Расчёт центров масс 
коленного сустава 


Алгоритм 
водораздела 


Вычисление точки 
каждого контура колена 


Вычисление 
и определение опорных 
линий и углов 


Рис. 6. Усовершенствованный алгоритм создания опорных линий и углов 
для обработки рентгеновских изображений коленного сустава 
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Обсуждение и заключение. Новый подход к исследованию рентгеновских снимков колена используется в 
больнице Аль-Базель (Латакия). Практика показала ценность усовершенствованного автоматизированного ана- 
лиза изображений для повышения точности диагностики. Отметим, что даже более ранняя версия программы 
эффективно выявляет аномалии кости, не замеченные при визуальном анализе [15]. С помощью предложенного 
решения ортопеды могут качественно прорабатывать малозаметные факторы, указывающие на нарушения, а 
также сочетания показателей, которые способны внести путаницу и помешать выявлению патологии. Из этого 
следует, что использование подхода, описанного в данной статье, позволяет обоснованно рассчитывать на повы- 
шение качества медицинских услуг в сфере ортопедии и хирургии. Специалисты получают инструментарий для 
более точного планирования и проведения операций. Кроме того, новое решение открывает возможности для 
лучшей организации консервативного лечения с учетом индивидуальных особенностей патологии. 

Следует продолжить научные изыскания в данном направлении для расширения диагностических возможно- 
стей и набора стратегий лечения заболеваний опорно-двигательного аппарата. Представляет интерес, в частно- 
сти, достижение абсолютной, т. е. 100-процентной точности алгоритмического определения линий. Следовало 
бы также сократить время обработки снимков. Оптимизация вычислительных процессов позволит получать ре- 
зультаты быстрее, что особенно актуально при неотложной диагностике. Кроме того, целесообразно доработать 
метод, чтобы применять его для анализа более сложных клинических случаев. Еще одно перспективное направ- 
ление — адаптация алгоритма для анализа других суставов, например тазобедренного или плечевого. 
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